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Abstract 


Flexible  masipulators  can  be  di2u:acterize<l  by  a  dynamic  model  with  a  large  number  of 
vibration  modes,  and  the  use  of  the  model  in  the  model-based  control  schemes  requires 
reduction  of  model  order.  Balanced  truncation  is  an  effective  method  for  model  reduction  of 
asymptotically  stable  systems  by  transforming  the  states  to  a  coordinate  system  in  which  the 
controllability  and  observability  Gramians  are  equal  and  diagonal,  and  eliminating  the  states 
which  contribute  weakly  to  the  input-output  map.  An  elastic  flexible  manipulator,  however, 
is  a  marginally  stable  system  and  thus  the  balanced  truncation  method  can  not  be  directly 
applied.  In  this  report,  we  present  a  method  of  reducing  the  order  of  a  marginally  stable 
system  bzised  on  the  fact  that  translation  transformations  in  the  frequency  domain  preserve 
input-output  properties  of  the  system.  We  address  the  successful  application  of  the  method 
to  model  reduction  of  flexible  manipulators  with  inflnite-dimensional  or  finite-dimensional 
model.  The  method  is  also  applicable  for  any  other  marginadly  stable  model,  such  as  elastic 
space  tnisswork  and  multi-dimensional  space  vehicle  structure. 
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1  Introduction 


Flexible  maniptilators  have  been  increasingly  interesting  to  many  researches  for  light-weight, 
energy-efficiency,  and  less  harmfulness  concerns.  Modeling  and  controlling  such  manipula¬ 
tors,  however,  is  a  challenge  because  conceptually  the  system  must  be  chairacterized  by  a 
distributed  parameter  systems  (DPS)  with  infinite  number  of  modes,  while  controlling  DPS 
requires  a  finite-dimensional  model.  A  various  approaches  of  methods  can  be  employed  to 
reduce  the  order  of  a  model  with  either  an  infinite-dimensional  or  a  finite-dimensional  struc¬ 
ture.  For  example,  model  order  can  be  reduced  by  retaining  the  first  finite  number  of  modes 
and  tnmcating  the  rest.  In  this  case,  the  criterion  for  selecting  appropriate  number  of  modes 
is  not  clearly  known  and  tnmcation  \isually  results  in  a  model  with  higher  or  lower  order 
than  that  is  needed.  Another  method  frequently  used  is  assuming  the  order  of  a  model  based 
on  lumping  mass  model  such  as  Lagrangian  approach.  Again,  this  method  also  has  difficiilty 
in  selecting  appropriate  mode  shapes  to  generate  an  accurate  model  with  a  minimal  number 
of  modes. 

The  balanced  tnmcation  method  [1]  is  based  on  the  realization  theory  that  for  zero  initial 
condition,  the  input-output  behavior  of  the  model  is  completely  specified  by  the  mutually 
controllable  and  observable  subspace  R'.  The  procedure  is  to  transform  state  space  descrip¬ 
tion  of  a  stable  system  to  balanced  coordinates  such  that  the  input-to-state  and  state-to- 
output  couplings  are  weighted  equally.  This  implies  that  observability  and  controllability 
Gramians  are  equal  and  diagonal.  The  diagonal  elements  of  the  balanced  Gramians  are  in 
fact  a  set  of  closed-loop  input-output  invariants,  or  the  Hankel  singular  values,  which  char¬ 
acterize  the  contribution  of  each  state  to  the  input-output  map  of  the  system.  The  states 
that  are  weakly  controlled  and  weakly  observed  have  ‘‘small”  singular  values  and  the  states 
that  contribute  weakly  to  input-output  map  can  be  ddeted. 

Model  reduction  by  balancing  techniques  is  applicable  only  for  stable  systems  [l]-[5].  This 
is  because,  for  the  balancing  techniques,  the  computation  of  the  controllability  and  observ¬ 
ability  Gramians  requires  that  the  system  be  asymptotically  stable.  Some  researches  [6]-[9] 
have  been  directed  to  applying  balancing  method  to  flexible  structure  by  assmning  lightly 
damped  and  widely  separated  modes.  This  assximption  is  normally  not  valid  in  practice. 
For  a  flexible  manipulator,  as  an  elastic  system,  we  have  to  face  the  fact  that  the  system  is 
marginally  stable.  There  are  very  few  research  articles  discussing  the  model  reduction  of  a 
marginally  stable  or  imstable  system.  In  the  paper  [10],  it  has  been  shown  that  an  unstable 
non-minimal  realization  (A,  B,  C)  can  be  transformed  via  a  similarity  transformation,  into 
a  balanced  one,  if  and  only  if  the  product  of  the  controllability  and  observability  Gramians 
is  similar  to  a  real  diagonal  matrix  A. 

This  report  presents  a  method  to  reduce  the  model  order  of  flexible  manipulators  using 
balancing  techniques.  The  approach  is  based  on  the  fact  that  translation  transformations 
in  the  frequency  domain  preserve  input-output  properties  of  the  system.  In  section  2,  we 
describe  the  modeling  of  two  flexible  manipxilators  with  em  infinite  or  finite  dimensional 
structure.  In  section  3,  we  present  the  theory  of  balancing  model  reduction  technique  for  a 
marginally  stable  system  where  the  poles  are  on  the  imajpnary  axis.  In  section  4,  we  address 
the  application  of  the  method  to  model  reduction  of  flexible  manipulators. 
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Figure  1:  The  One  Link  Flexible  manipulator 

2  Modelling  of  Flexible  Manipulators 


In  this  section,  we  will  discuss  the  modeling  of  flexible  manipulators  with  infinite-dimensional 
and  finite-dimensional  structures.  These  two  models  will  be  used  to  demonstrate  the  proce¬ 
dure  of  model  reduction  in  section  4. 


2.1  Infinite-dimensional  model 


The  infinite-dimensional  model  for  a  one-link  flexible  manipulator  shown  in  Figure  1  can 
be  modeled  [11]  under  the  following  assumptions  :  (1)  motion  is  planar  and  payload  and 
joints  are  symmetric  objects  in  XY  plane.  (2)  rotary  inertia  and  shear  deformation  effects 
are  neglected.  (3)  gravitation  force  is  in  the  Z  direction.  (4)  beam  inertia  and  flexibility  are 
uniformly  distributed  over  the  link  length. 

Based  on  these  assumptions,  the  dynamics  of  the  position  of  any  point  along  the  link  is 
governed  by  the  Bemoulli-Euler  beam  equation: 


d*y{x,t)  .  ^y{x,t) 


=  0 


(1) 
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where  0  <  x  <  I  and  the  equation  has  associated  boundary  conditions. 

The  transfer  function  for  general  BemouUi-Eulet  beam  equation  can  be  found  in  [12]  and 
the  complete  dynamics  can  be  described  by 


- +“—^^ =/(*.<) 


dt^ 


Q(x,0)  =  Qo(x),  ^(x,0)  =  Q,(x),  |^(0,i)  =  in(t) 

0(0.0  =  «.(<),  0(<.O  =  S3(i),  0(/,«)  =  j4(t) 

0  <  X  >  I,  t  >  0,  a  >  0 


(2) 


To  solve  the  above  problem,  we  may  define  a  standardizing  function: 

tn(a:,<)  =  /(x,t)  +  (:?o(x)^'(t)  +  -  aV(a:)^i(t)  +  o^5(a:)if2(<)  “  (3) 

-aH{l  -  x)g3(t)  -  a^S(l  -  x)g4(t) 

and  we  convert  the  original  problem  into  a  zero  initial  and  boundary  condition  problem  for 
which  the  Green’s  function  is 


<Pn(x)(pn(0 


sin  an\t 


(4) 


where 


¥’n(x)  =  {sink  fi„l  —  sin  finl){cosh  finX  +  cos  HnX) — 
—  {cosh  Uni  —  cos  flnl){sinh  flnX  +  StTl  /t„x) 

and  the  fin  are  the  non-negative  roots  of  the  equation  cosh  filcos  fil  =  1. 
The  transfer  function  is  given  by: 


which  has  infinite  poles: 


«  <^n(0  +  OVJI 


s  =  ±jafil,  n  =  l,2,... 


(5) 


2.2  Finite-dimensional  model 

For  a  multi-degree  of  freedom  flexible  manipulator,  Lagrangian  dynamics  is  an  appropriate 
approach  to  model  the  system.  As  an  example,  the  Self- Mobile  Space  Manipulator  developed 
in  the  Robotics  Institute  of  Carnegie  Mellon  University  [13][14]  can  be  modeled  as  shown  in 
Figure  2.  Two  flexible  links  are  connected  by  the  flex  joint  and  the  end  of  the  link  holds  the 
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Parameters 

li 

i»l-3 

length  of  link  ■' 

m; 

is  2. 3 

mass  of  joint  i 

J. 

i>l*3 

joint  inertia  i  with  high  gear  ratio,  (/,»» 

A|*:fcarr«Uo 

/„^-;«MUrliNnla 

bending  stiffness  of  link  i,  (X’^-  s  —Lj.' 
Cl  t’ 

MCU— wiabwrtU 

torsional  stifhwss  on  link  i.  (JCi,-  =  — 7- 

/g :  palir  SMMMM  cT  w  iMTtto 

Kbi 

i«l3 

Ka 

•  »13 

Variables 

e< 

is  J-3 

angtilar  disf^emeni  of  joint  i 

V,i,  Va 

i-2.3 

bending  defleokm  at  the  end  of  link  1 

i-2.3 

torsional  deflection  of  link  i 

♦ji-  ♦zi 

i»2.3 

angular  bending  deflection  of  link  i 

i=I-3 

lonjue  on  joint! 

Figure  2:  The  Two  Link  Flexible  Manipulator 
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tip  mass,  while  the  other  end  of  the  link  is  connected  to  the  flex  and  twist  joints  attached 
to  the  base.  Each  compact  joint  contains  a  driving  motor  and  a  high  reduction  gearing,  and 
each  flexible  link  is  slender  tubing  to  maximize  its  strength  with  minimum  weight. 


The  following  assumptions  are  made  to  simplify  the  model:(l)  The  middle  joint  and  the 
end-effector  are  lumped  as  a  point  mass  and  the  moment  of  inertia  of  the  motor  in  the  joint 
is  considered  for  the  high  reduction  gearing.  (2)  The  mass  and  the  moment  of  inertia  of  the 
links  are  neglected,  compared  to  the  joints  and  tip  mass,  or  the  moment  of  inertia  of  the 
motor.  (3)  Each  link  possesses  high  tensile  and  shear  stiffness  so  that  no  tensile  and  share 
distortion  are  assumed.  Compliance  in  the  bending  and  torsional  directions  of  the  link  is 
considered.  (4)  The  deflections  of  the  link  is  small  with  respect  to  the  length  of  the  link.  (5) 
The  gravity  effect  is  not  taken  into  account  for  space  applications. 

The  kinetic  energy  is  a  summation  of  the  rotationad  kinetic  energy  of  three  revolute  joints 
and  the  translational  kinetic  energy  of  two  point  mass.  The  potential  energy  consists  of 
torsional  strain  energy  and  bending  strain  energy.  Based  on  the  kinetic  energy  and  potential 
energy,  Lagrangian  model  is  as  follows. 


where 


i=2 


-?L  ?IL-  f. 

dt  dqi  dqi  ^  dqi  ’ 


q  =  [  •  •  •  gi3  j 

—  [  ^2  VjQ  ^V2  ^*3  ^3  ^113  ^*3  ] 

f  =  [  /l  •  •  *  /l3  ] 

=  [  2i  Ta  0 . 0  T3  0 . 0 


(6) 


(7) 

(8) 


The  joint  angular  position  and  link  deformation  in  torsional  and  bending  directions  are 
selected  as  a  set  of  generalized  displacements,  and  the  joint  torques  and  zero  elements  form 
a  set  of  generalized  forces.  By  a  standaird  procedure  based  on  Lagrangian  dynamics,  the 
explicit  dyn2imics  equation  can  be  obtained. 


M{q)q  +  K  (q)q  +  C{q,  q)  =  f 


(9) 
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where  M  is  the  inertia  matrix,  K  is  the  stiffness  matrix,  and  C  is  the  centrifugal  and  Coriolis 
torques. 

The  dynamic  equation  above  is  nonlinear  due  to  the  translational  kinetic  energy  part  in 
Equation  (7).  To  linearize  the  model  with  reasonable  accuracy,  we  assume  that  the  speed 
of  the  robot  is  low  so  that  centrifugal  and  Coriolis  torques  are  negligible  and  the  effect  of 
link  deformation  on  the  variation  of  the  inertia  matrix  is  small  so  that  the  inertia  matrix  is 
a  function  of  only  joint  angles.  In  this  way  the  d3mamic  model  can  be  linearized. 

Mq  +  Kq  =  /  (10) 

It  is  obvious  that  the  poles  of  the  model  lie  on  the  ju  axis. 

The  model  can  be  decoupled  into  two  independent  parts,  tangential  model  and  radial  model. 
The  tangential  model  represents  the  motion  associated  with  joint  1  driving  the  tip  mass  in 
the  circular  direction,  while  the  radial  model  describes  the  motion  associated  with  joints 
2  and  3  driving  the  tip  mass  in  the  plane  perpendicular  to  the  circular  direction.  It  must 
be  noted  that  the  mass  and  stiffness  matrices  for  both  models  are  functions  of  the  mixed 
joint  angles.  In  this  way,  we  can  obtain  two  linearized  models  in  the  tangectial  and  radial 
directions 

Mtqt'  -I-  Ktqt  =  Tt  (11) 


Mt 


Kt 


ji  +  cos^  ^2  +  ^  cos  $2  +  4  4 
m£  cos  ^2  +  4  m  +  1  1 

4  1  1 


qt  =  [ 

=  I  ^  i  (u*3  +  4  sin  +  £3  cos  03^^ )  ]  ^ 

Tt  =  [  fi  0  0 
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/x  =  +  cos((l2  +  fla)  fc  =  ki  = 

A>s^  +  =^  +  J.  "  =  \/5  t'-t 

Mrqr"  +  =  Tr 

j2  +  +  ^  +  1  +  2£  cos  ^3  m£  +  ^  +  cos  63  1  +  £  cos  63  1  +  £  cos  63 


ni£  +  £  4"  cos  ^3 

m  + 1 

cos  03 

cos  ^3 

1  +  £  cos  63 

cos  63 

1  + is 

1 

1  +  £  cos  ^3 

cos  ^3 

1 

1 

0  0  0  0  • 

0  ijfc<y+3)  Q 

0  0  0  0 

0  -iShi  0 

L  ^2  A2  -i 

[  ^2  h  UV3 

[<?2  ^  <?3 

[  f2  0  73  0  ]^ 

[^  0  ^ 

32  =  h  =  ;^, 

A2  s  4jfc£2  +  3 


(12) 
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3  Model  Reduction  for  a  Marginally  Stable  System 

FVom  the  discussion  above,  it  has  been  known  that  the  model  of  flexible  manipulator  is 
m2urginally  stable,  i.e.,  the  poles  of  the  system  lie  on  imaginary  axis.  In  this  section,  we 
present  a  method  of  model  reduction  for  such  a  system. 

Consider  the  linear  time  invariant  system 

x(t)  =  Ax{t)  +  Bu{t) 

y{t)  =  Cx{t)  (13) 

where  x(t)  €  il"  is  the  state  vector,  u(t)  €  iT”  is  the  input  or  the  control  vector,  i/(t)  6  BF 
is  the  output  vector,  and  A,  B,  C  are  matrices  of  appropriate  sizes. 

The  transfer  function  of  the  system  (1)  is  given  by 

H{s)  =  C{sI-A)-^B  (14) 

Definition  1:  The  symbol  {A,B,C)  represents  a  state  space  realization  of  the  transfer 
function  of  the  system. 

Definition  2:  The  realization  (A,  B,  C)  is  said  to  be  minimstl  if  it  is  completely  controllable 
and  observable. 

Definition  3;  {A,B,  C)  is  said  to  be  asymptotically  stable  if  Re  (  A  |  A  |)  <  0. 

Definition  4:  Let  (A,  B,  C)  be  a  minimal  and  asymptotically  stable  realization.  (A,  B,  C) 
is  said  to  be  baljmced  if  A,  B  and  C  satisfy  the  matrix  Lyapunov  equations  with  equal  and 
diagonal  controllability  and  observability  Gramians,  that  is 

AWc  +  WcA^  =  -BBi^  (15) 

A^Wo  +  WoA  =  -C^C,  (16) 

where  Wc  =  Wo  =  ^  =  diag[(ri  •  o-„}  and  Hankel  singular  values  <Ti  >  (Tj  >  03  >  . . .  >  0 


The  matrix  S  can  be  partitioned  into  two  submatrices  Ei  and  E2  in  the  following  way 


E,  0 

0  Ej 


(17) 


where  Ei  =  diag[<Ti  aj  •  ■  •  a*]  and  Ej  =  diag[ak^\  <Tjt+2  •  •  •  <Tn]- 


Further  partition  A,  B,  C  corresponding  to  the  partition  of  E 


A  = 


An  Ai2 
A21  A22 


(18) 
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Then  the  kth  order  reduced  order  model  is  Ci). 

Definition  5:  A  minimal  realization  {A,B,C)  is  said  to  be  a  J  form  realization  if  A  is  a 
Jordan  matrix  J  and  ( J,  B,  C)  has  following  structure: 

J  =  diag[Jx  J2  •  •  •  Ja]i  B  =  [Bn  Bn  *  *  *  Bia]^ 

C  =  [CnCn-‘-Ci^]  (19) 

with 


and 


Jij  —  diag[Jij^  Jij^  Bij  —  Bij^  Bij^^.\, 

Cij  =  [Ci j,  Ct ] 


Jijk  — 


Ai  1  0  •••  0  0 

0  A.  1  0  •  •  •  0 

0  0  •*.  **.  0 

:  :  **.  A.  1  0 

:  :  •••  0  Ai  1 

0  0  .  0  Ai 


dij\  —  «4j2  —  *  •  •  =  JijTijt  Jijk  €  C 
Bijk  —  [^tjffei^ijfcj  •  •  Bijk  € 

^ijk  —  [A'iAriCijfcj  ’"^jkai}]i  ^ijk  €  C  . 


If  (A,  B,  C)  is  a  minimal  realization,  there  exists  a  non-singuletr  matrix  T  such  that 


J  =  T-^ATy  B  =  T-^B,  C  =  CT 


(20) 


H{s)=C{sI-J)-^B  (21) 

In  the  rest  of  this  section,  we  present  an  approach  to  model  reduction  of  marginally  stable 
system  using  balancing  technique.  We  at  first  state  the  following  theorems. 

Theorem  1:  The  controllability  and  observability  of  a  minimal  realization  J-form  {J,  B,  C) 
depends  only  on  the  structures  of  J,  B  and  C. 

The  proof  of  Theorem  1  can  be  obtained  directly  from  the  definitions  of  the  controllability 
and  observability. 

Theorem  2:  For  minimal  marginally  stable  realization  H{s),  there  exist  translation  trans¬ 
formations  in  frequency  domain,  s'  =  s  -1-  a,  such  that  the  poles  of  H{s  -|-  a)  lie  strictly  on  the 
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left  half  s-plajie  and  the  controllability  and  observability  of  H{s  +  a)  axe  identical  to  those 
of  H{s). 

Proof:  If  H(s)  is  a  minimal  realization,  we  have 

H{s)  =  C(sI-J)-^B 


and 


H(t)  = 

Let  Xmax  be  the  largest  eigenvalue  of  the  matrix  J.  For  the  marginally  stable  system 
R^{Xmax)  =  0.  we  can  choose 

a  =  e 

where  e  >  0. 

A  A 

Let  H{t)  =  then  the  Laplace  transform  of  H{t)  can  be  determined  by 

His)  =  His  +  a)  =  C[(s  +  a)I  -  =  C(s/  -  J)-^B 


Thus 


ReiX  I  J  I)  <  0 


A  A  A 

Since  rank(  J)  =  rank(  J),  ( J,  B,  C)  is  a  minimal  realization  and  J  and  J  have  same  structure. 
Prom  Theorem  1,  ( J,  B,  C)  and  ( J,  B,  C)  have  the  same  controllable  and  observable  subspace. 


In  the  balanced  realization,  the  ith  Hankel  singular  value  characterizes  how  controllable 
and  observable  the  ith  state  is.  Therefore  the  model  can  be  reduced  by  neglecting  the 
states  corresponding  to  “small”  Hankel  singular  values.  This  implies  that  the  criterion  for 
balanced  model  reduction  depends  only  on  the  input-output  information  of  the  system.  Thus 
it  is  possible  to  obtain  a  reduced  order  model  by  balancing  method  for  an  asymptotically 
stable  system  if  the  input-output  map  of  the  system  can  be  determined.  More  specifically, 
we  transform  the  marginally  stable  system  into  the  asymptotically  stable  system  by  the 
translation  transformation  in  the  frequency  domain  and  balance  the  resulting  stable  system. 
This  procedure  leads  to  the  following  algorithm. 

Algorithm 


A  A 

1.  Transform  His)  into  His)  =  ff(s  +  a)  so  that  His)  has  no  poles  on  the  right  half 
s-plane  nor  on  the  ju  axis. 

2.  Find  a  minimal  state  space  realization  (yl,  B,  C)  of  His). 

3.  Compute  the  balanced  realization  iAiyBi„Ci,)  from  iA,B,C).  That  is 

Ai  =  T-'^AT,  Bi^T-^B,  Ci,  =  CT, 

where  T  is  the  transformation  matrix  which  makes  the  controllability  and  observability 
Gramians  equal  and  diagonal. 

T  can  be  obtained  using  any  of  methods  presented  in  [l]-[5]. 


10 


4.  Determine  the  reduced  order  model  (A,,  Br,  Cr)  from  (Aj,,  Bt,  Cb).  Recall  that  for  the 
balanced  realization,  W'c  =  IVo  =  S  =  diag[ai  oj  •  •  •  On]  and  Hankel  singular  values 
O’!  >  0-2  >  o’a  ^  >  0*  If  o’fc  then  S  can  be  pcjtitioned  as 

V  _  f  Si  O' 

^  [  0  Ss 

where  Si  =  diag[<cTi  <72  •  •  •  Ok]  and  £2  =  dtaiif[o‘it+i  o-*+2  •  •  •  o‘n]- 


(Aft,  Bft,  Cb)  can  be  partitioned  accordingly 


Aa  = 

Bu  = 

'  Bb,  ■ 

)  jjb 

,  a  =  [Cft,  Cft,] 


to  yield  the  reduced  order  model 

Af  —  Aftj^  ,  Bf  —  Bbi ,  Cr  ~  • 

5.  Obtain  Hris)  from  (Ar,J5r,Cr)  as 

Hris)  =  CrisI-Ar)-^Br. 

6.  The  reduced  order  model  of  the  original  system  is  given  by 

Hris)  =  H{s  -  a). 


Some  remarks  may  be  in  order: 

1.  The  translation  transformation,  s'  =  s  +  a,  is  a  conformal  transformation  which  pre¬ 
serves  congruence  and  involves  neither  rotation  nor  stretching.  Therefore  the  translation 
transformation  in  the  frequency  domain  does  not  affect  the  input-output  description  of  the 
model. 

2.  The  algorithm  can  be  interpreted  in  the  time  domain  as  follows 

Hit)  =  Hit)  •  e-“*  ♦  e“‘  =  Hit)  •  e“‘  (22) 

Hrit)  =  Hrit)  •  e-“* .  e»‘  =  Hrit) .  e“‘,  (23) 

A 

where  Hit)  and  Hit)  are  the  impulse  responses  of  the  full  order  model  and  the  modified 
full  order  model,  Hrit)  and  Hrit)  are  the  impulse  responses  of  the  reduced  order  model  and 
the  modified  reduced  order  model,  respectively.  Elquation  (23)  suggests  that  the  impulse 
response  of  an  unstable  system  Hit)  can  be  factored  into  a  stable  impulse  Hit)  multiplied 
by  c“*,  while  Equation  (24)  implies  that  Hrit)  is  equal  to  Hrit)  multiplied  by  c“‘.  Therefore, 
if 

Hrit) «  Hit) 

then 

Hrit) «  Hit). 
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3.  Some  model  reduction  methods,  such  as  Routh  approximation,  may  be  applicable  to  the 
unstable  model  reduction,  but  they  may  result  in  a  stable  reduced  order  model  and  malce 
implementation  of  model-based  control  scheme  difficult.  For  the  proposed  method,  the  poles 
which  lie  on  the  imaginary  axis  will  be  the  dominant  poles  of  the  modified  system  and  thus 
balancing  will  keep  properties  of  these  poles. 

4.  The  method  can  be  applied  equally  well  to  both  SISO  and  MIMO  systems.  In  MIMO 
case,  H{s)  is  a  transfer  function  matrix. 
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4  Model  Reduction  of  Flexible  manipulators 


As  discussed  in  section  2,  the  flexible  naanipulator  can  be  modeled  as  either  inflnite-dimensional 
model  or  finite-dimensional  system.  In  either  case,  it  is  generally  needed  to  reduce  the  order 
of  the  model  so  that  any  model-based  control  scheme  can  be  implemented  in  practice.  In 
this  section,  we  discuss  the  application  of  the  method  proposed  in  the  previous  section  to 
the  two  models  developed  in  section  2. 

Example  1:  Infinite-dimensional  model 

We  list  five  modes  of  the  infinite-dimensional  model  based  on  the  parameters  suggested  in 
[11]  in  Table  1. 


Table  1:  Poles  and  zeros  of  the  flexible  manipulator 


Pair 

Poles 

Zeros 

1 

0.00 

±11.18 

2 

±15.42; 

±60.45 

3 

±49.97; 

±149.27 

4 

±104.25; 

±277.58 

5 

±178.27; 

±445.35 

This  is  a  tenth  order  model  and  is  difficult  to  be  implemented  in  real-time  control.  Using 
our  method,  we  can  reduce  the  model  order  and  the  resultant  second  order  model  is 

5^  +  1.6013  X  lO^s  -  7.6801  x  10® 

~  s2  -h  3.3671  X  10-3s  -  3.3552  x  lO-* 

The  step  responses  of  the  tenth  order  model  and  the  corresponding  second  order  model  are 
compared  in  Figure  3.  It  has  been  shown  that  the  second  order  model  is  nearly  identical  to 
the  tenth  order  model. 

Example  2:  Finite-dimensional  model 

The  second  example  of  flexible  manipulators  is  the  radial  direction  model  developed  in  section 
2  which  has  two  inputs  and  four  outputs  and  is  manipulator  configuration  dependent.  At  62 
=  03  —  60  degree,  the  transfer  function  is  given  by 

NUMn  =  3*(2.6823e  -I-  Ols^  +  3.7052c  +  03s^  -I-  9.1806c  -|-  03) 

NUM12  =  s’(1.7131c  +  03s2  -  1.3589c  +  04) 

NUM2i  =  s^(-2.6823c  -I-  Ols^  -  3.5937c  -I-  03) 

NUMt2  =  -1.9382c  -I-  03s^ 

NUM31  =  a’(1.7131c  +  035’  -  1.3589c  +  04) 

NUM32  =  s*(7.4510c  +  Ols^  -I-  7.6274c  -1-  033^  +  3.0999c  +  04) 

NUMn  =  s^(-2.6823c  -I-  Ols^  -  5.4970c  +  03) 
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—  Original  model  . . .  Reduced  order  model 


Figure  3:  Step  Responses  of  the  Full  and  Reduced  Order  Models 
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NUMu  =  5^(-7.4510e  +  Ols^  -  9.1002e  +  03) 

DEN  =  s\s*  +  2.3010c  +  025^  +  I.1032e  +  04) 

where  NUMij/DEN  represents  the  transfer  function  from  the  torque  t,-  to  the  angle  $j.  This 
is  a  8th  order  MIMO  system  and  using  our  method  the  reduced  to  7th  order  is  ^l5  follows. 

NRii  =  3.4686e  -  01s«  +  2.6394e  +  01s®  +  7.6040e  +  Ols^  +  3.6210e  +  03s®  +  3.3186e  +  03s®  + 
7.0192c  +  03s  -  9.9127c  +  02 

NRi2  =  1.6638c  -  01s«  -  3.4956e  -  01s®  +  3.8582e  +  Ols**  +  1.6517c  +  03s®  +  1.6239e  +  03s®  - 
1.4633c  +  04s  +  2.0666c  +  03) 

NR21  =  s®(-4.6450c  -  04s^  -  2.6773c  +  01s®  +  3.7313s®  -  3.5851D  +  03s  +  5.0513c  +  02) 

NR22  =  s®(-2.2280c  -  04s*  +  2.4323c  -  02s®  -  6.9621c  -  02s®  -  1.9340c  +  03s  +  2.7257c  +  02) 

NR31  =  1.6667c  -  01s«  -  2.5671c  -  01s®  +  3.8477c  +  Ols^  +  1.6660e  +  03s®  +  1.6091c  +  03s®  - 
1.4628c  +  04s2.0660c  +  03) 

NR32  =  7.9945c  -  02s®  +  7.4317c  +  01s®  +  7.9596s^  +  7.5946c  +  03s®  -  1.7942c  +  02s®  + 
3.0496c  +  04s  -  4.3073c  +  03) 

NR^i  =  s®(-4.1424c  -  04s^  -  2.6699c  +  01s®  +  3.6426s®  -  5.4776c  +  03s  +  7.7092c  +  02) 

NR42  =  s®(-1.9869c  -  04s‘‘  -  7.4450c  +  01s®  +  1.0504c  +  01s®  -  9.0908c  +  03s  +  1.2826e  +  03) 

D£JVie  =  a®(s®-1.4289c-01s^+2.3010c+02s®-3.2765c+01s®+1.1032c+04s-1.5582c+03) 
where  NRijfDENR  represents  the  transfer  function  of  the  reduced  order  model  from  the 
torque  r,-  to  the  angle  6j.  The  step  responses  of  the  original  model  and  the  reduced  order 
model  are  compared  in  Figure  4. 
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—  Original  model  . . .  Reduced  order  model 


Figure  4:  Step  Responses  of  the  Full  and  Reduced  Order  Models 
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5  Summaries 


We  presented  a  simple,  eflBcient  method  for  model  reduction  of  flexible  manipulators  based  on 
the  fact  that  the  translation  transformation  in  s-plane  preserves  input-output  property  of  the 
system.  By  translation  transformation  in  frequency  domain,  we  can  change  the  stability  of 
the  system  vdthout  lose  of  its  input-output  property.  The  original  b2dancing  method  requires 
that  the  model  be  asymptotically  stable  while  this  method  is  applicable  to  a  marginally  stable 
system.  We  employed  the  proposed  method  to  reduce  the  order  of  flexible  manipulator 
models  such  that  model-based  scheme  can  be  implemented.  The  method  has  been  examined 
by  two  examples,  an  inflnite-dimensional  model  and  a  finite-dimensional  model. 
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